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CALCULATION OF AERODYNAMIC FORCES ON AN 
INCLINED DUAL-ROTATING PROPELLER 
By John L. Crigler and Jean Gilman, Jr. 


SUMMARY 


A method is herein presented for calculating the aerodynamic forces 
on the blades of a dual-rotating propeller which has its thrust axis 
inclined to the air stream. This method makes use of Theodorsen's pro- 
peller theory for the calculation of the aerodynamic forces on the pro- 
peller blade. At any particular blade station and blade position the 
instantaneous flow field is presumed to be a steady-state condition for 
the whole propeller. Application of the method involves the use of two- 
dimensional airfoil characteristics as determined in a steady flow. 

Calculations were made for an 8-blade NACA 8.75-(5) (05)-037 dual- 
rotating propeller operating at advance ratios of 1.20 and 5-25 for 
flight Mach numbers of 0.30 and 0.623, respectively. The thrust axis 
angle is 4° in both cases. Results of the calculations show that, in 
general, the fluctuations in aerodynamic forces on the front component 
are greater than the fluctuations in forces on the rear component. Com- 
parisons with the forces on a single -rotating propeller show that the 
magnitude of the fluctuating forces on the front compone- of a dual- 
rotating propeller is of the same order as the forces occurring on a 
single propeller operating at the same advance ratio, blade-angle setting, 
and inclination of the thrust axis. 

Inasmuch as the front component of a dual-rotating propeller under- 
goes about the same variation in thrust coefficients as would a single- 
rotating propeller, the fact that smaller variations are encountered on 
the rear blades means that the one-cycle-per-revolution vibration problem 
for a pitched or yawed dual-rotating propeller is no more serious than for 
the single -rotating propeller. 


ffc- 


2 


NACA RM L53D30 


INTRODUCTION 


In reference 1 the magnitude of the aerodynamic forces on a single- 
rotating propeller in pitch or yaw was calculated under the assumption 
that the existing propeller theory may he used as though steady- state 
conditions existed successively at several a.ngul a.r positions of the pro- 
peller blade during one revolution. The method made use of the usual 
steady-state compressible -airfoil characteristics with Goldstein correc- 
tion factors for a finite number of blades. The calculated results were 
compared with experimental data taken from a pitched propeller installa- 
tion. The magnitude of the aerodynamic forces on the propeller blade 
calculated frcsn the steady-state method was found to be in good agree- 
ment with the experimental data although there was some lag in the forces 
as shown by the experimental data. Further calculations based on the 
theory of osci lla ting airfoils in a pulsating flow (ref. 2) showed that 
the lag in air forces found in the experiments could be predicted. 

The procedure followed for the case of the single- rotating propeller 
is inadequate in the dual case. In addition to the unsteady flow in the 
propeller plane caused by pitch or yaw, the effect of the flow field of 
each component on the other component for each position of the blade must 
be considered. In the case of the dual- rotating propeller the flow field 
of each component has been considered in making calculations of the forces 
on both components of a dual -rotating propeller in an effort to predict 
the magnitude of the differences in forces on the front and rear compo- 
nents. In view of the good agreement obtained in the case of the s ingle - 
rotating propeller, it was assumed that steady- state conditions existed 
successively at the angular positions of the dual -rotating propeller 
during each revolution. 

The oscillating air forces on the propeller blades cause vibratory 
stresses having a frequency of one per revolution. In this report the 
air forces are determined for each component of the dual-rotating pro- 
peller in an effort to determine the relative severity of the one -cycle - 
per-revolution vibration problem on the two components and to make a 
comparison with one component of the dual-rotating propeller when 
operating as a single propeller. 


SYMBOLS 


B number of blades 

b 
c d 


chord of propeller-blade element, ft 
section drag coefficient 
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dCrp 

dx 

h 

J 

K(x) 

M 

n 

R 

r 

T 

t 

V 

W 

Wo 

v 

W 

X 

a 

C&p 

p 


section lift coefficient 
design lift coefficient 

thrust coefficient, t / 
propeller diameter, ft 


element thrust coefficient. 


dT/dx 
pn^D^ - 

blade- sect ion maximum thickness, ft 
advance ratio, V/nD 
circulation function 
Mach number of advance 
propeller rotational speed, rps 
tip radius, ft 

radius to any blade element, ft 


thrust, lb 
time, sec 

velocity of advance, ft/sec 

resultant velocity at blade section, ft/sec 

geometric resultant velocity, ft/sec 

rearward displacement velocity of helical vortex sheet 

velocity ratio w/V cos 

radial location of blade element, r/R 

angle of attack, deg 

angle of inclination of propeller thrust axis, deg 
blade angle, deg 
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K 


tan"*^- 


£d 

C I 

mass-flow coefficient. 


PI 


JO 


K(x)x dx 


p mass density of air, slugs/cu ft 

a solidity, Bb/2itr 

ac-i propeller-element load coefficient 

(f> aerodynamic helix angle, deg 

0 O geometric helix angle, deg 

a> angular velocity of propeller, 2jm, radians /sec 

Subscripts : 

F front 

R rear 

S. front component of a dual-rotating propeller operating as a 

single-rotating propeller 

cot denotes angular position of blade required to give instantaneous 

values, deg 


AERODYNAMICS OF THE INCLINED FROFELLER 
The Velocity Diagram 


The velocity diagram for a blade section of a single-rotating 
inclined propeller is discussed in reference 1. This diagram is repro- 
duced herein as figure 1. The convention adopted in reference 1 was to 
consider the time variable cot, which defines the angular position of 
the blade, to be zero when the blade Is initially vertically upward with 
the propeller axis In positive pitch arid measured in the direction of 
rotation (the axes used In defining the blade position may be rotated 
to comply with propeller attitudes other than pitch). 
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From figure 1, the geometric helix angle is given "by 


(H* ' ta “' 1 


V cos cty 

imDx + V sin ctq sin cut 


or 


Wo* nx 


cos 


= tan" 


OCip 


-j- + sin Orp sin ojt 


The resultant velocity is given by 


(W°)ctft = \ffi cos^ctrp + ^itnDx + V sin oq sin cut) ‘ 


The local advance ratio is given by 


( 1 ) 


( 2 ) 


J cut 


«x cos oq 


JtX 

J 


+ sin otrp sin cut 


( 5 ) 


From, equation (3) the local advance ratio is seen to vary, the amount of 
variation depending on the position of the blade. In making steady-state 
calculations of the forces and moments on the blade of a single-rotating 
pitched propeller in reference 1, a change in blade position was treated 
as a change in operating V/nD in accordance with equation (3)- For the 
single -rotating propeller the solution is adequate when the calculations 
are based on the local geometric helix angle inasmuch as complete calcu- 
lations were made as though the entire propeller operated successively 
at the different blade positions and the induced effects were incorporated 
in each calculation. 

The velocity diagram for the blade sections of a dual-rotating pro- 
peller showing the effect of the induced flow of each component on the 
other component is given in figure 2. This diagram is reproduced from 
reference 3« 

For inclined dual-rotating propellers the alteration to the geometric 
velocity field due to oscillating velocities for the front propeller unit 
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is similar to the alteration for single-rotating propellers, that is, 


tan 



V cos cftp 

JthDx + V sin c&p sin cot 


± ( 4 ) 

22. + sin Cbj\ sin cot 

The rear propeller is of opposite rotation and., in order to make use of 
the same reference axes for measuring the blade-position angle as is 
used for the front propeller, it is convenient to write 


tan (<MH 


V COS Orj, 

JthDx - V sin arp sin cot 


cos ctij 


JtX 


sin dip sin cot 

J X 


( 5 ) 


This reversal of sign for the front- and rear-propeller units is used 
throughout the development of the equations for the aerodynamic helix 
angles and element loading coefficients. 

In addition to the alteration to the geometric velocity field for 
the dual-rotating propeller, it is necessary to consider the influence 
of the velocity field of each component of the propeller on the other 
component before performance calculations can be made. Theodor sen’ s 
propeller theory was used to determine the flow field of each component 
due to the other component. Admittedly, this method has the defect that 
it was developed for the idealized case in which the two components have 
the same load, distribution and equal power absorption. However, if it 
is assumed that the mass-flow coefficient and the circulation functions 
for a given number of blades apply for nonoptimum distribution, the method 
with these assumptions gives a close approximation to the flow field 
since both of these functions are primarily dependent on the geometric 
helix angle and show only small variations with loading for a given 
geometric helix angle and blade number. This dual-rotating -propeller 
case is analogous to the use of Goldstein functions for the several blade 
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positions in the case of single -rotating propellers as developed in 
reference 1. The interference velocities for the dual-rotating propeller 
suid the effect of each component on the flow field of the other component 
are shown in reference 5* Since the front and rear components of the 
dual-rotating propeller are unequally loaded for each phase angle and 
the loading changes with a change in phase angle, it is necessary to 
make successive approximations to the mass-flow coefficient k and to 
the displacement velocity w for each phase angle calculated. 


Method of Calculation 

The resultant velocity, the aerodynamic helix angle, and the element 
loading coefficient for a radial station of an optimum dual-rotating 
propeller are given in reference 3 as 


= idrS 1 + * KV sln2 ^ 

Wr ■ TirrS 1 + * " elaS *°) 


(6) 

(7) 


tan rftr = — 
c nx 


1 + i w^l + i k tan 2 0^ 


tan | w^l - | K tan 2 0) 



irx 


(l + w)W sin 0 o K(x) 
1 + i kv sin 2 0 o 



(l + v)w sin 0 o K(x) 
1 + ^ kw sin 2 0 o 


C8) 

(9) 

( 10 ) 


itx 


( 11 ) 
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where 


ton *- si 1 + 1 5 ) (12) 

A method of applying Theodorsen's propeller theory to the calcula- 
tion of the aerodynamic forces on the propeller "blade for nonideal dual- 
rotating propellers is described in reference 4. The method essentially 
consists of solving for local values of w along the blade by a method 
of cross-plotting by using equations ( 8 ) to ( 12 ). 

The method, as developed in reference 4, has been found to be satis- 
factory when the loading on the two components of the propeller is close 
to the same but should not be used when the difference in loading becomes 
large. Consequently, the c hang es developed in the present report have been 
incorporated in the addenda to reference 4. The method, as developed 
herein, shows good agreement with experimental data even for large dif- 
ferences in loadings on the two components of the propeller as shown in 
the addenda of reference 4. 

For inclined dual -rotating propellers the propeller-performance equa- 
tions require modification to account for the oscillating velocities. The 
resultant velocity, the aerodynamic helix angle, and the element loading 
as modified are given to a first approximation for the front components by 


V COB On / i _ 0 

( W Ojt)p = + 15 ; K 03t W 03t sin rC 


'cut 


(13) 




, V COS CCrn 

, = tan -1 

imDx + V sin Qip sin cut 


( 1 *) 


COS 


*“(*«*). ■ 




1 + 1 "tot! 1 + 1 v tan 2 0 ] 


(15) 


— + sin cu p sin cat 
J A 
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where 


W 03t = 


V CSt 

V cos oti p 


tan 0^ 


V COS Q£rj-i( 1 + 



jtnDx + V sin c&p sin cut 




(16) 


The effect of the inclined thrust axis on a dual -rotating propeller is 
seen to be the same as rotating the two components at unequal speeds. 
Equations for the rear component are obtained by a similar modification. 


In calculating propeller characteristics by Theodorsen's theory as 
applied in reference 3 the average axial interference velocity kw was 

assumed to be equally due to each of the two components. When the 

loading is not the same on the front and rear components, or effectively, 
the rotational speeds are unequal, as in the present case, the inter- 
ference velocity will not be equally due to each component and may be 

looked upon as broken into two unequal parts (kw) f and (kw) r where 

the average axial interference velocity in the fi n al wake is 


KW 


(kw) f + (kw) r 
2 


but 


(kw) f 0 (kw) r . The expression (kw) f refers 


to the interference velocity that would be obtained in the final wake 
if the loading on the rear component were equal to the loading on the 
front component of a dual-rotating propeller. Similarly, (kw) r refers 

to the interference velocity that would be obtained if tl„ loading on 
the front component were equal to the loading on the rear component of 
a dual-rotating propeller. Thus, when the two components are equally 
loaded and of opposite rotation, the interference velocity is equally 
due to each component and kw = (kw) f = (kw) r . 


In the solution the interference velocity associated with each 
component will result in different values of 0 ^ for the two components 

of the propeller. The values of w F and w R are first calculated 

separately just as though each component operates in the presence of 
another component equally and oppositely loaded. The method is equiv- 
alent to assuming that ideal functions apply to nonoptimum loadings and 
ideal values of k and K(x) are used. This practice has been generally 
successful In single-rotating propellers where Goldstein's tip correc- 
tion factors are used for nonoptimum loadings. In the solution it is 
adequate to replace 0 by 0o F or 0^ in the Induced velocities. 

Equations (13) and ( 15 ) as modified now become 


10 


NACA KM L 55 DJO 


V COS C&P 

TJ . L J 

S~ 

1+1 

(^cutp + 

Bin 2 f i] 1 

F 

1 + S 

2 



( 17 ) 





1 ^ 


— w 


a&F 


+ i(icv) 


Ojtp 



— + sin atj\ sin cut 


(18) 


where 





JtX 

J 


+ sin o&p sin cat 


( 19 ) 


Similarly, for the rear propeller with due regard to the reversal in 
sign because of opposite rotation in obtaining the aerodynamic helix 
angle, 


W(0b R 


V COS OCrp 

B ln ( 0 o ) ajfcR 



^)fl5t g + C^MtH 
2 




i (in '>a<t p B±n2 (fc) aJ t E 

or 


V COS CSrp 

^ 8in (A>)eBtR 



( 20 ) 
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and 


tan ^cut R = 


cos 


+ \ v cut R " |(* w )cut F tan 2 ^! + JUw)^ tan 2 ^^ 


— - sin ay sin cut 
J 


( 21 ) 


where 


tan tfVutp = 


cos 


+ | *CUt R ) 
2 itx . . 

— - sin Off sin cut 


(22) 


It can he seen that, for equal loadings on the front and rear compo- 
nents where (tcw) F = (kw) r , equations ( 17 ), (l8), (20), and (21) reduce 

to the same form as in reference 3- 

The element loading coefficient as modified for the inclined dual- 
rotating propeller is 



cos ay 

(1 + sin(0 o ) aJ t/(x) ajtF 


/ nx . 

1 — + sin ay 

sin cut^ 

i + ^| 

( KV CUt F 

V 2 I 

sin 2 (0 o )ajt F 


( 23 ) 


where [kU)^^, is based on - = *x tan^o)^^! + v^) for the 

front component and 


( ao iW 


cob 0^1 + v (0tR jw 0JtR aiK^JojtpKfx)^ 


_ sin ^ sin ojfcUi + | sin 2 0 O ^ 


r( KV W + r( KW ) 


cut R 


( 24 ) 

where (^(x)^ R is based on icx tan (0o) ^tR (l + w cut R ) for the 
component. Finally, the thrust equation is given by 


rear 
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dT = cos 0 - c^ sin 0 )dr 


( 25 ) 


By definition. 


x 


_ £ _ 2r 
R D 


dx = 


2dr 

"E" 


a = £b_ 

2rtr 

c, 

tan 7 = — 
c l 

therefore, equation ( 25 ) may be reduced to 

2 

— = -^2— xpW^cc, cos $f( 1 - tan 7 tan 0) (26) 

dx 4 

By substituting the value of W from equation (17) and replacing 0 
by the value of 0^ and dividing both sides by pn 2 D^ the thrust 
equation for the front component in coefficient form is reduced to 


idx / , 

\ / ojtp 


XX jgc ° aga T 

T sln^) 


ost F 






^2 


Bin 2 (^ 0 ) 


OjtW 


cos( ^ajt F E ' tBn(9>) art F ttm ^mt F 


( 27 ) 
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Similarly, for the rear component 


W /ojfcR 


nx ^cos 2 ^ r 

T sln2 ( 0 °U E L 


+ 


1 

2 


sin2 (^o) 


Oit R 





(“iW c ° s((!> >o* E [^ 


- tan^)^ tan 



( 28 ) 


RESULTS OF SAMPLE CALCULATIONS 


Calculations of the elemental- thrust variations have been made for 
an 8 . 75 -foot -diameter 8 -blade dual- rotating propeller having an 
NACA 8.75-(5) (05)-037 blade design. Blade-form curves for this propeller 
are given in figure 3* NACA l 6 -series airfoil data appropriate to the 
various section Mach numbers encountered in making the calculations were 
interpolated from reference 5 * 

In order to Include the effect of advance ratio on the fluctuating- 
forces calculations were made for both a low advance ratio (j = 1 . 2 ) and 
for a relatively high advance ratio (j = 3*25) • The thrust axis angle 
Orji was maintained constant at 4°. This angle was selected as being 
large enough to show significant effects of Inclination on the fluctuating 
forces but not so large as to create blade-section angles of attack In 
the stalled range where airfoil data are usually unavailable for making 
the calculations. For the low-speed case (j = 1.2) the flight Mach 
number is 0.30. The blade -angle settings are fJp Q ^ = 32.41° and 

Prq ^ = 31 . 53 °. For the high-speed case (J = 3 * 25)7 the flight Mach 

number is 0.623* The blade-angle settings are Pp-Q ^ = 57*89° and 

(3rq = 56 . 15 °* la either case the blade -angle settings used would 

give equal power absorption on the front and rear components with = 0°. 

Figures 4 and 5 show typical results for several radial stations inves- 
tigated for J = 1.2 and J = 3*25, respectively. Instantaneous values 

of Cm . as a function of cut are shown In figures 6 and 7 for 
■ L 0 jtF j R 

J = 1.2 and J = 3*257 respectively. The elemental-thrust coefficients 
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for the front component of the dual-rotating propeller operating as a 
single -rotating propeller at the same J, cuj>, and £q ^ were calcu- 
lated by the method given in reference 1 for one condition in order to 
give a comparison of the difference in forces on a single- and dual- 
rotating propeller. These results (for J = 5.25) are included in fig- 
ure 5 WTld the integrated values showing the variation of the instantaneous 
thrust coefficient Crp^ 3X6 shown in figure 7- 

Both the elemental-thrust coefficients for the various radial sta- 
tions plotted as a function of angular blade position (figs. 4 and 5) 
and the integrated thrust (figs. 6 and 7) show that the difference between 
the maximum and minimum forces, shown by the fluctuations in thrust coef- 
ficients, is greater on the front component than on the rear component 
of a dual propeller. The fluctuations in thrust coefficients on the 
front component of a dual propeller, due to pitched operation, differ 
only slightly from the fluctuations in thrust coefficients on a single- 
rotating propeller operating at the same J, ctrp, and blade-angle setting. 

The value of (Ct^^qo - ^=2700^ f ° r th * lov " B P eed caBe 

fig. 6) is 0.0420 for the front propeller and is 13 percent higher than 
the value of “ ^^=90°) R ° f °‘°^ for the rear propeller. 

A similar comparison for the high-speed case (J = 3-25> fig- 7) giveB a 
value of (^^.cjqO " G T a jfc=270°) f of which is 28 percent higher 

than the value of (Cm „ - Cm . __ Q \ of 0.168 for the rear pro- 

\ i ojt=270° 1 ojt=90°/R 

peller. The comparison of the fluctuation in forces of the front 
component of the dual-rotating propeller with the single -rotating pro- 
peller would be expected to show about the same variations since the 
velocity field of the front component is affected by the average axial- 
interference velocity only of the rear component which is, in general, 
rather small. The rear component is affected by the axial- and rotational- 
interference velocities of the front component, but the rotational- 
interference velocity tends to decrease the magnitude of the fluctuating 
forces on the rear component compared with the front component. Inasmuch 
as the rotational interference velocities increase with an increase in 
operating J, the fluctuating forces on the rear component decrease 
compared with the forces on the front component. 

This discussion is only intended to apply to the oscillating air 
forces on the blades that cause vibratory stresses having a frequency 
of one per revolution. Other effects such as blade passage and effects 
of blade structure are not within the scope of this paper. 
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SUMMARY OF RESULTS 


A method is herein presented for calculating the fluctuating aero- 
dynamic forces on pitched or yawed dual-rotating propellers. Typical 
examples have indicated the following results: 

1. Comparison of the fluctuating aerodynamic forces due to pitched 
or yawed operation of a dual-rotating propeller show that, in general, 
the fluctuation in forces on the front component are greater than the 
fluctuation in forces on the rear component. 

2. The differences in the magnitude of the fluctuating forces 
increase as the operating advance ratio is increased. 

J. The magnitude of the fluctuating forces on the front component 
of a dual-rotating propeller are of the same order as the forces 
occurring on a single -rotating propeller operating at the same advance 
ratio, blade-angle setting, and inclination of the thrust axis. 

4-. The one-eycle-per-revolution vibration problem for the dual- 
rotating propeller is no more serious than for the single -rotating pro- 
peller when operating at the same advance ratio, blade-angle setting, 
and inclination of the thrust axis. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va. 
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Design lift coefficient, 





U>t, 


(c) X = 0.9« 

, /dC T \ 

Figure 4.- Typical results of calculations shoving I i as a 

Vdx /afi F)R 

function of cut for several radii of the NACA 8*75- (5) (05)-037 
propeller. M = O.JOj J = 1.2; 0 &p - 4°. 
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(b) x = 0.7. 

( ac T \ 

Figure 5 .- Typical results of calculations showing I- — as 

V 4 * Mf,R 

function of cut for several radii on the NACA 8.75-(5) ( 05 ) — 03T 
propeller. M = 0 . 623 ; J = 3 « 25 ; ccm - 4°. 
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if instantaneous thrust coefficient 
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